Quasi-linear transport approach to equilibration of quark-gluon plasmas 
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We derive the transport equations of quark-gluon plasma in the quasi-linear approximation. The 
equations are either of the Balescu-Lenard or Fokker-Planck form. The plasma's dynamics is as- 
sumed to be governed by longitudinal chromoelectric fields. The isotropic plasma, which is stable, 
and the two-stream system, which is unstable, are considered in detail. A process of equilibration 
is briefly discussed in both cases. The peaks of the two-stream distribution are shown to rapidly 
dissolve in time. 



I. INTRODUCTION 

The quark-gluon plasma (QGP), which is produced in 
relativistic heavy-ion collisions, is believed to be equili- 
brated within a time interval of order of 1 fm/c or even 
shorter [lj. Such a fast equilibration is naturally ex- 
plained assuming that the quark-gluon plasma is strongly 
coupled 0-13. Then, scattering processes are very fre- 
quent and relaxation times are short. However, the the- 
ory of high-energy density QCD [j| suggests that due to 
the existence of a large momentum scale Q s , at which 
the gluon density saturates, the plasma is rather weakly 
coupled at the early stage of the collision because of 
asymptotic freedom. Experimental data on jet quench- 
ing indicate that the coupling constant a s < 0.3 6, 7], 
even though the value assumes averaging over the whole 
evolution of the QCD medium created in the relativis- 
tic heavy-ion collision. Thus, the question arises how 
fast the weakly interacting plasma equilibrates. Due 
to anisotropic momentum distributions the early stage 
plasma is unstable with respect to the chromomagnetic 
plasma modes. The instabilities isotropize the system 
and thus speed up the process of its equilibration. The 
scenario of the instabilities-driven isotropization is re- 
viewed in [8]. However, the complete evolution of the 
plasma momentum distribution is now accessible only by 
numerical simulations [9l4llj|. 

The transport theory of a weakly coupled quark-gluon 
plasma has been studied since 1980s when the kinetic 
equations in the mean-field approximation were derived 
[HCDl. Although the mean-field dynamics is rather sim- 
plified, the equations are still difficult to solve due to their 
non-linear structure. If one is interested in small devi- 
ations from equilibrium or any other homogeneous and 
stationary state, the equations can be linearized and then 
solved. The mean- field transport theory, which is lin- 
earized in small deviations from equilibrium, is now well 
understood, for a review see [14j . It is known to be equiv- 
alent to the effective QCD in the hard-thermal loop ap- 
proximation. The linearized transport theory around any 
homogeneous and stationary but non-equilibrium plasma 



state was also worked out and the connection with the 
dia gram matic hard loop approximation was established 
Tol 16| . Numerous problems of the theory of the quark- 



gluon plasma were successfully resolved within the hard 
loop approach. For example, a systematic method to 
eliminate infrared divergences, which plague perturba- 
tive calculations, was developed, see the reviews [l7l.[l8|. 

However, various questions cannot be addressed within 
the mean-field theory. For example, transport coefficients 
are then formally infinite. Thus, there were numerous ef- 
forts to derive transport equations of quark-gluon plasma 
which hold beyond the hard loop approximation |19l - [39l | . 
These efforts were mostly concerned with the transport 
properties of an equilibrium quark-gluon plasma. Our 
motivation is rather different. We are interested in equi- 
libration of quark-gluon plasma, in particular in the equi- 
libration of the system which is initially unstable. Thus, 
we intend to study how fluctuating deviations from a 
quasi-stationary non-equilibrium state influence the sys- 
tem's bulk or average momentum distribution. This ef- 
fect of back reaction is particularly important in the case 
of unstable systems. The linear response theory describes 
how unstable modes initially grow in the presence of 
a non-equilibrium momentum distribution, but it says 
nothing on how the modes modify the plasma momentum 
distribution. Thus, the problem of equilibration cannot 
be addressed in such a theory. 

Our objective here is to derive the transport equations 
where the bulk distribution function slowly evolves due 
to the interaction with fluctuating chromodynamic fields. 
We actually consider only a simplified problem of QGP in 
a self-consistently generated longitudinal chromoelectric 
field. This simplification is not much needed for isotropic 
plasma but it appears crucial to study anisotropic sys- 
tems. Taking into account only the longitudinal chromo- 
electric field, we obtain the transport equations of the 
Fokker-Planck or Balescu-Lenard form which describe 
the effect of back reaction. A similar, but incomplete 
effort was undertaken by Akkelin (39|. The derivation 
presented here closely follows the procedure developed 
for the electromagnetic plasma, where it is known as the 
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quasi-linear theory or the theory of a weakly turbulent 
plasma [45-42] . The theory assumes that the distribu- 
tion function of plasma particles can be decomposed into 
a large but slowly varying regular part and a small fluc- 
tuating or turbulent one which oscillates fast. The av- 
erage over the statistical ensemble of the turbulent part 
is assumed to vanish and thus the average of the distri- 
bution function equals its regular part. The turbulent 
contribution to the distribution function obeys the colli- 
sionless transport equation while the transport equation 
of the regular part, which is of our main interest here, 
is determined by the fluctuation spectra. The fluctua- 
tions of chromodynamic fields, which are used to derive 
the quasi-linear transport equations, were studied in [43| 
where stable and unstable plasma states were considered. 

The Fokker-Planck equation derived here is somewhat 
similar to the equation obtained in [44], EH . It was used 
there to show that the chromomagnetized quark-gluon 
plasma exhibits an anomalous shear viscosity, as pres- 
ence of the domains of chromomagnetic field leads to the 
momentum transport in the plasma. 

Our paper is organized as follows. In Sec.|TI]we present 
the QGP transport equations; the notation and conven- 
tions are introduced. The decomposition of the distri- 
bution functions into the regular and turbulent parts is 
discussed in Sec. IIIII The explicit expressions of the 
fluctuating distribution functions which obey collision- 
less transport equations are derived in Sec. lIVl A general 
form of the equations of the regular distribution functions 
is found here as well. Further discussion splits into two 
parallel parts: Sec. [V] is devoted to the stable isotropic 
plasma while in Sec. IVII the unstable two-stream system 
is discussed. Although we neglect transverse chromody- 
namic fields, the collision terms of transport equations, 
which are found here for an isotropic plasma, are very 
similar to those derived in [H HH, ^ Hi, [H, . As 

an application of the transport equations we derived, a 
process of equilibration of the isotropic plasma and of the 
two-stream system is discussed. The paper closes with a 
summary of our considerations and outlook. 



II. PRELIMINARIES 

The transport theory of a quark-gluon plasma, which 
forms the basis of our analysis, is formulated in terms of 
particles and classical fields. The particles - quarks, anti- 
quarks and gluons - should be understood as sufficiently 
hard quasiparticle excitations of quantum fields of QCD 
while the classical fields are highly populated soft gluonic 
modes. An excitation is called "hard" , when its momen- 
tum in the equilibrium rest frame is of order of the tem- 
perature T, and it is called "soft" when the momentum 
is of order gT with g being the coupling constant. Since 
we consider a weakly coupled quark-gluon plasma, the 
coupling constant is assumed to be small g <C 1. In our 
further considerations the quasiparticles are treated as 
classical particles obeying Boltzmann statistics but the 



effect of quantum statistics can be easily taken into ac- 
count. 

The transport equations of quarks, antiquarks and glu- 
ons are assumed to be of the form 

DQ(i,r,p)-i{F(f,r),V p Q(t,r,p)} = 0, 

T>Q(t,r,p) + i{F(i,r),V p Q(*,r,p)} = 0, (I) 

0G(t,r,p)-i{F(t,r),V p G(t,r,p)} = 0. 

The (anti-)quark distribution functions Q(t, r,p) and 
Q(t,r, p), which are N c x N c hermitean matrices, be- 
long to the fundamental representation of the SU(7V C ) 
group, while the gluon distribution function G(t, r, p), 
which is a (JVf — 1) X (JV| — I) matrix, belongs to the ad- 
joint representation. The distribution functions depend 
on the time (t), position (r) and momentum (p) vari- 
ables. There is no explicit dependence on the time-like 
(/i = 0) component of the four-vector p^ as the distri- 
bution functions are assumed to be non-zero only for 
momenta obeying the mass-shell constraint p^p^ = 0. 
Because the partons are assumed to be massless, the ve- 
locity v equals p/E p with E p — |p|. T> = D° + v • D is 
the covariant substantive derivative given by the covari- 
ant derivative which in the four-vector notation reads 
£)/.( = Qfj. _ ig^A^^x), • • • ] with A^(x) being the chromo- 
dynamic potential. The mean-filed terms of the transport 
equations ([T]) are expressed through the color Lorentz 
force F(t,r) = g(~E(t,r) + v x B(<,r)). The chromo- 
electric E(i, r) and chromomagnetic B(i, r) fields belong 
to either the fundamental or adjoint representation. To 
simplify the notation we use the same symbols V, D°, 
D, E, and B for a given quantity in the fundamental 
or adjoint representation. The symbol {...,...} denotes 
the anticommutator. 

The collision terms are neglected in the transport equa- 
tions fT}. The collisionless equations are applicable in 
three physically different situations: when the distribu- 
tion function is of (local) equilibrium form; when the 
timescale of processes of interest is much shorter than 
the average temporal separation of parton collisions; and 
when the system dynamics is dominated by the mean 
field. In our study we refer to all three situations. When 
the equilibration of isotropic plasma is discussed, it is cru- 
cial that the collision terms vanish in local equilibrium. 
In the case of unstable two-stream system, the effects 
of collisions can be initially neglected, as the growth of 
unstable modes is very fast. Later on, the strong fields 
become mostly responsible for the system's evolution 

The transport equations are supplemented by the 
Yang-Mills equations describing a self-consistent gener- 
ation of the chromoelectric and chromomagnetic fields. 
The equations read 

D-E(t,r) = p(t,r) 
D • B(i, r) = , 

D x E(*,r) = —D B(t, r) , (2) 
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D x B(t,r) - j(i,r)+A)E(t,r) , 

where the color four-current j 7 ' = in the adjoint 

representation equals 



^^lr 



(2.)3 Sp - L ^,r,p) 
T Q (Q(;,r,p)-Q(t,r,p)) 



(3) 



where r a , T a with a = 1 , . . . , iV c 2 - 1 are the SU(iV c ) group 
generators in the fundamental and adjoint representa- 
tions, normalized as Tr[r a r b ] = \5 ab and Tr[T Q T b ] = 
N c 6 ab . The set of transport equations (QJ and Yang- 
Mills equations ([2]) is covariant with respect to SU(iV c ) 
gauge transformations. 



III. REGULAR AND FLUCTUATING 
QUANTITIES 

We assume that the chromodynamic fields and distri- 
bution functions which enter the set of transport equa- 
tions can be decomposed into the regular and fluctuat- 
ing components. The quark distribution function is thus 
written down as 



Q{t,r,p) = (Q(t,r,p))+8Q(t,r,p) 



(4) 



where (• • •) denotes ensemble average; (Q(t, r, p)) is 
called the regular part while SQ(t, r, p) is called the fluc- 
tuating or turbulent one. It directly follows from Eq. (|4]) 
that (SQ) = 0. The regular contribution is assumed to 
be white, and it is expressed as 



<Q(t, r,p)> =n(t,r,p)I 



(5) 



where I is the unit matrix in color space. Since the distri- 
bution function transforms under gauge transformations 
as Q — > U QU^ 1 , where U is the transformation matrix, 
the regular contribution of the form ([5]) is gauge inde- 
pendent. We also assume that 



1(0)1 » \SQ\ 

but at the same time 



|V P (Q)| » \V P SQ\ , (6) 



\VSQ\ > \V(Q)\ . (7) 



Analogous conditions are assumed for the antiquark and 
gluon distribution functions. What concerns the chromo- 
dynamic fields, we assume in accordance with ([5]) that 
their regular parts vanish and thus 



85Q 




d(Q) 


dt 


> 


dt 



<E(t,r)) = (B(t,r)>=0. 



(8) 



We substitute the distribution functions (Q} into the 
transport equations and the Yang-Mills equations and 
linearize the equations in the fluctuating contributions. 



The linearized transport and Yang-Mills equations re- 
main rather complex. Therefore, we discuss here a sim- 
plified problem: we consider a QGP in the presence of 
turbulent longitudinal chromoelectric fields, but neglect 
the chromomagnetic and transverse chromoelectric fields. 
This simplification can be avoided for an isotropic plasma 
but it is needed, as explained below, to make progress on 
an analytical treatment for anisotropic systems which are 
our main interest here. The simplified transport equa- 
tions then read 

Z><fQ(t,r,p)-ffE(t,r).V p n(t,r,p) = 0, 
V5Q{t,v,p) + gV{t,v) -V p n(t,v,p) = 0, (9) 
2>aG(t,r,p)-0E(t,r)-V p n fl (t,r,p) = 0, 

where T> = + v • V denotes from now on the material 
(not covariant) derivative. 

The equation describing the self-consistent generation 
of a longitudinal chromoelectric field is 

V • E„(t, r) = p a (t, r) = -g J -0^ SN a {t, r, p) , (10) 

where 

<W„(t,r,p) = Tr[T a (5Q(t,r,p)-SQ(t,r,p)) 

+T a SG(t,r,p)] . (II) 

The linearized equations are formally Abelian but they 
include a fundamentally non-Abelian effect, i. e. the 
gluon contribution to the color current. Therefore, the 
gluon-gluon coupling is partly taken into account. The 
linearized Yang-Mills equation corresponds to the multi- 
component electrodynamics of -Y c charges (in the so- 
called Heaviside-Lorentz system of units). The equations, 
however, are no longer manifestly covariant with respect 
to SU(iV c ) gauge transformations. Nevertheless, our final 
results are gauge independent. 

We now substitute the distribution functions (U]) into 
the transport equations ([1]). Instead of linearizing the 
equations in the fluctuating contributions, we take the 
ensemble average of the resulting equations and trace 
over the color indices. Thus we get 



Vn-^-Tr(E-\7 p 5Q) 
2?n+^Tr(E- V P SQ) 
Vn g -j^-jTr(E-V p SG) 



(12) 



Since the regular part of distribution function is assumed 
to be color neutral, see Eq. ([5]), the terms of the form 
Tr[(E • V p n)] vanish because the field E is traceless. The 
trace over color indices also cancels the terms originating 
from covariant derivatives like Tr([A M , SQ]). We finally 
note that the regular distribution function n is gauge 
independent and so is Tr(E • V p SQ). 
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IV. SOLUTION OF THE LINEARIZED 
EQUATIONS 

Due to the condition (0, the space-time dependence 
of the regular distribution functions is neglected in the 
linearized transport equations © and the equations be- 
come easily solvable. We solve Eq. (j9|) with the initial 
conditions 



JQ(t = 0,r,p) - <5Q (r,p) , 
*Q(t = 0,r,p) = <5Q (r, P ) , 
<5G(i = 0,r,p) - <5G (r,p) , 



(13) 



using the one-sided Fourier transformation denned as 

/(w,k) = y dt J d 3 re 4(w *- k r) /(i,r) . (14) 

The inverse transformation is 
D+la duj f d 3 k 



-oo+ia 27r J (27T) 



e- i ^ t - k ' r V(w,k) , (15) 



where the real parameter a > is chosen in such a way 
that the integral over uj is taken along a straight line in 



the complex lu— plane, parallel to the real axis, above all 
singularities of /(cj,k). We note that 



iw/(w,k) = /(t = 0,k) 



(16) 



The linearized transport equations ©, which arc 
transformed by means of the one-sided Fourier transfor- 
mation, are solved as 



6<2(w,k,p) = i 

lu — n. • v 

srit u \ . ffE- V p n(p) -(5Q (k, p) 
dQ(u>,k,p) = i — — — — , (17) 

<5G(w,k, p) = i 



gE ■ 


V p n(p) + 


^Qo(k,p) 




lu — k 


V 


gE ■ 


Vpfl(p) - 


SQ {k, p) 




lu — k 


V 


gE ■ 


V p n 3 (p) - 


MG (k,p) 




cj k 


• V 



We note that the color electric field E(w, k) retains its full 
frequency and wave number dependence in these equa- 
tions. Inverting the one-sided Fourier transformation, 
one finds the solutions of linearized transport equations 
as 



6Q(t,r,p) = g[ dt' E(t', r - v(f - 0) • V P n(p) + SQ (r - vt,p) , 
Jo 

5Q(t, r, p) = -g f dt' E(t', r - v(i - 0) ' V P n(p) + SQ (r - vt, p) , (18) 
Jo 

6G(t,r,j>) = gf dt'E(t',r-v(t-t'))-V p n g (p) + 5G (r-vt,p), 
Jo 

where we assumed that E(lu, k) is an analytic function of lu. With the help of solutions (p~8|) . the force terms in the 
transport equations (|12[) become 



(E(t,r)-V p( 5Q(t,r,p)) = g [ dt' V p (E\t, r)E* (f, r - v(i - 0)>Vjn(p) + V P r)<5Q (r - vt, p)) , 

Jo 

(E(*,r) .V p 5Q(t, r,p)> = -g [ dt' %{E\t, r)E* (f , r - v(t - 0)>Vjn(p) + V l p (E l (t, r)SQ (v - vt, p)> , (19) 

Jo 

(E(t,r).Vp«G(i,r,p)> = <? / di' V p (£ l (t, r)E* (f, r - v(i - f)))V p n 9 (p) + %{E\t, r)SG (r - vt, p)> , 

Jo 



We conclude that the transport equations (IT^|) are determined by the correlation functions (E' l (t ) r)E 3 ] (t\ r')), 
r)(5Qo(r', p)), (-E 4 (t, r)<5Qo( r ', p)), an d r)<5Go(r', p)). To compute these functions, the state of the plasma 

must be specified. Although we are mainly interested in an anisotropic plasma, we start with the isotropic case. 
Thereafter, we consider the two-stream system. 
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V. ISOTROPIC PLASMA 



For the case of isotropic plasma, the correlation functions of both longitudinal and transverse fields are well known 
[43j . Here we limit our considerations to longitudinal chromoelectric fields, whose correlation function is [43j : 

2 f-ca+ia 1 r- 73. /.oo+ier' j / c- j3 y 

(EUt,r)Ei(t',r')) = 9 -8<* / ^ " ^M+.V-^-k'-O (20) 



2 y_ 00+iCT 2^ j (2tt)3 y_ oo+iff , 2^ y (2^ 

k*k' 3 (27r) 3 (5( 3 )(k' + k) /■ d 3 p /(p) 



x 



k 2 k' 2 ei(w, k) e £ (w', k') 7 (2tt) 3 (w - k ■ v)(w' - k' • v) ' 



where /(p) = ra(p) + n(p) + 27V c n £ ,(p) and Sl(uj,Il) is 
the longitudinal chromodielectric function discussed in 
the Appendix. Note that we do not assume that n(p), 
n(p), and n g (p) are given by the thermal equilibrium 
distributions, only that they are isotropic functions of p. 
Zeroes of sl(uj, k) and £l(w', k') as well as those of the 
denominators (lo — k • v) and (&/ — k' • v) contribute to 
the integrals over u> and u) 1 . However, when the plasma 
system under consideration is stable with respect to lon- 



gitudinal modes, all zeroes of el lie in the lower half-plane 
of complex u). Consequently, the contributions associated 
with these zeroes exponentially decay in time, and they 
vanish in the long time limit of both t and t' . 

We are further interested in the long-time limit of 
{E l (t, r)E l (t', r')). The only non-vanishing contribution 
corresponds to the poles at to = k • v and u' = k' • v. 
This contribution reads 



{E a {t,v)E b {t , r )) - 2 d y (27r)3 y (27r) 3 e k 4 | £i ( k .v,k)p ' (21) 

The correlation functions like (E l a (t, r)SQo(r' , p')) are not computed in ref. [43[ but they can be readily inferred 
from the formulas given there. One finds 

(EHt r)6Q (r> p')> - -«r° f"*' *L f **- f ^-^-^ *' (^) 3 ^ (3) + k) n( P Q 

- laT - f d " k e -ik.(v't-r +r Q) fc' "(pO f22l 

" 5 7 (27T) 36 k 2 e L (k-v',k) ' 

where the last equality holds in the long-time limit which is carried by the contribution corresponding to the pole 
uj = k • v'. Similarly, one finds 

^ du f d 3 k f <Pk> „_ iM _ k . r _ k ,. fc' (2^ 3 )(k' + k) n(p>) 



{El(t,r)SQ (r' lP ')) = gr a / — / — . / -— 



+ia 2irJ (2tt) 3 7 (2tt) 3 k 2 £ L (w,k) w - k • v' 

" ff y k 2 e L (k-v',k) ' (23) 

(EHt r)SG (r> r>')) aT« [°°*" d0J f f ^ e -(^'^' (2Tt)H^{U + k) n g (p') 

- jnT a ( p -ik-(v't-r+r')) n g{p') 

Substituting the correlation functions (f2"Tl |2"2"|) into (TlT)|) . one finds 

T,.(E( ( , ,) ■ W r, P)> - £tf - 1) jf *V, / ^ / <f|e— »"'>^ Vj»(p) 
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and analogous expressions for Tr(E(t, r) • V p 5Q(t, r, p)) and Tr(E(i, r) • V p SG(t, r, p)). As shown in |43j . 
Tr(E' l (t, r)E 3 (£', r')) is gauge independent within the linear response approach. The same arguments used to show 
this apply to Tv(E t (t,r)SQo(r' ,p')). Thus, we conclude that the collision term of the transport equation (| 1 2 [1 . 
Tr(E(£,r) • V p SQ(t, r, p)), is gauge independent. 

Let us now discuss the first term on the r.h.s. of (|25j) . Computing the integral over t' we get 



Tr<E(i,r).V P .5Q(*,r,p)) (1) = ^(iV c 2 - 1) V* J -0^ J 



d 3 k k l k j 



/(p') 



x — i 



(2tt) 3 J (2tt) 3 k 4 \e L (k- v',k)| 2 
cos(k • (v - v')t) - 1 sin(k • (v - v')t) 



(26) 



k • (v - v') 



k- (v 



V p n(p) 



The first term does not contribute to the integral because it is an odd function of k. Since in the limit t — > oo we 
have 



sin(k • (v — v')t) . , , v , 

lim ) / -ri-L = k • v - v' 

t^oo k ■ (v - v') v ' 



one finally finds 



Tr<E(i, r) • V p 5Q(t, r, P )) (1) = £ n (N 2 C - 1) V£ J 0L J ^ 



d 3 k k l ti 



Z(p') 



i 3 k 4 |e L (k-v',k)| 



(27) 



5(k-(v-v')) V>(p).(28) 



Analogously, one computes Tr(E(i,r) • \7 p SQ(t, r, p))(i) and Tr(E(i,r) • \7 p 5G(t, r, p))(i) • 
The second term on the r.h.s. of (1251) can be written as 



Tr(E(t, r) • V p SQ(t, r, p)) (2) = §(A 2 - l)Vj, 



d 3 k k l 3e L (k-v,k) 



(2tt) 3 k 2 |e L (k-v,k)| 2 



n(p) 



(29) 



because the term with 5?£z,(k • v, k) is an odd function of k (see Appendix). Alternatively, one can argue that 
the r.h.s. of (|2"!)|) has to be real as the l.h.s. is real. In the same way one finds: (E(t,r) • V p 6Q(t, r, p))(2) and 
(E(f,r)-V p <SG(t,r,p)) (a) . 

With the formulas derived above, the transport equations (|12|) can now be written either in the Balescu-Lenard 
form or the Fokker-Planck form. ^ 



A. Balescu-Lenard equations 



Using the formula (|A.2[) to express Sel through the 
distribution function, the transport equations (|12[) get 
the Balescu-Lenard form [42j 

V n(t, r,p) = V p • S[n, n,n g ] , 
T>n(t, r, p) = V p • S[n, n, n g ] , (30) 
Vn g (t,r, p) = V p • S g [n,n, n ff ] , 

where, as previously, T> is the material derivative, and 
d 3 p' 



5" 



7 '<<) "-gj 



S l [n,n,n g 



Vjn(p)/(p')-n(p) Vj,/(p') 

^^(v V) 
(2tt) 3 1 ' J 

V^(p)/(p')-n(p)V^/(p') 

(2tt) 3 ff 1 ' J 
V p n g (p)/(p')-n g (p)V^/(p') 



(31) 



with 

y , _ ff 4 A 2 - 1 /• rf 3 fc fc'fcJ 27r<5(k • (v - v')) 

(V ' V) - 8 A c J (2tt) 3 k 4 |e L (k-v,k)| 2 ' 

(32) 

and 



B?(vy) 



2A e 2 
A 2 -l 



(33) 



Since the interaction processes that are taken into ac- 
count conserve the numbers of particles of every species 
(^iQiS); t ne transport equations in the Balescu-Lenard 
form (|30[) can be seen as continuity equations in momen- 
tum space with S, S, S g playing a role of currents. One 
observes that for classical equilibrium functions 

-E p /T 



/ cq (p),n c >),7l c >)X q (p) ~ 



(34) 

the collision terms (|3"Tj) vanish, as expected, because 

( 1) , -/)B y (v,v')=0. (35) 

If £i(o;,k) is replaced by unity, i.e. if one ignores 
the chromodielectric properties of the plasma, the ten- 
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sor S lJ (v, v') is easily found to be 

oil, i\ ^ N ^- 1 L 



32tt N c |v-v'| 



(36) 



x S* 



with 



L 



(v - v') 2 



dk/k = ln(fc max /fc 

mm J • 



(37) 



The parameter L is called the Coulomb logarithm and 
the collision term with the tensor B lJ '(v 1 _v') of the form 
(f36|) is called the Landau collision term |42j. Estimating 
fcmax as the system temperature T and fc m ; n as the Debye 
mass itid ~ one finds L ~ ln(l/<?). 

It may appear strange that we start with the collision- 
less transport equations (JlJ to derive the collision terms. 
This procedure, which is commonly used in the plasma 
literature, is well justified, however, see e.g. Ref. [42j . 
The collision terms, which are derived above, represent 



J 



the effect of fluctuating soft fields on the hard quasipar- 
ticles. It is important to note that the collision terms are 
dominated, as it should be, by the soft wave vectors. 
Consequently, the collisions of quasiparticles involving 
the exchange of hard momenta, which are neglected in 
Eqs. (p}, do not need to be taken into account at lowest 
order. 



B. Fokker-Planck equations 

Sometimes it is more convenient to use the trans por t 
equations in the Fokker-Planck form. Following Ref. [42| , 
one rewrites Eqs. (|3U)l as 

(P-V^WVj-V^Wjnftr.p) = 0, 
(2>-Vp^>)Vj-VjY i (v))n(t,r,p) - 0,(38) 
(V - V;*«(v)V£ - V;y;(v))n fl (t, r, p) = , 

where 



^ c 2 -i) 



d 3 P ' r d 3 k fc<# Z(p') . p/, , ,^ 

d 3 k k l 3e L (k-v,k) 
(^k 2 " |e L (k-v,k)P 

dV f d 3 k V k-V p ,/(p') 



|^/( P ')^(v,v'), (39) 



(40) 



d 3 p' 



(27T) 



V£,/(p')i^(v,v') 



3 P 



and 



iV? - 1 



X«(v), 

y'(v) . 



(41) 
(42) 



The equations (|38j) appear to be linear but actually they 
are not: the coefficients X ij (v), F l (v), X«(v) and Y*(v) 
depend on the distribution functions. When the dis- 
tribution functions are of the classical equilibrium form 
(/ eq (p), ™ eq (p), n eq (p), n<p(p) ~ e~ E ^l T ), we have the 
relation 



y<(v) = -x«(v) . 



(43) 



Consequently, the Fokker-Planck collision terms vanish 
in equilibrium, as do the Balescu-Lenard collision terms. 

Since the system is assumed to be isotropic, X 1 ^ (v) 
and Y l (w) can be expressed as follows: 



X ij (v) 



a6 l 



(44) 



nv) 



with 



b = - 



d 3 p' 
(2tt) 3 

(2tt) 3 



/(pO 
Z(p') 



3vV l - <F 1 



(45) 

5 u '(v,v') (46) 
B y '(v,v') (47) 



« i Vi/(p')S tf (v,v') 



(48) 



Because of the system's isotropy, the coefficients a, 6, 
c can depend only on v 2 . In the ultrarelativistic limit, 
which is adopted here, v 2 = 1, and consequently a, b, c 
are independent of v. We also note that in equilibrium 
the coefficients are related as 



T c = a + b , 
which follows from Eq. (|43|) . 



(49) 
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When £i(o;,k) is replaced, as previously, by unity one 
finds that 6 = and 



96tt 



■(N 2 c -l)L / dpp 2 f(p) 
Jo 



96tt 



<N 2 c -l)L 



dpp' 



df(p) 
dp 



(50) 



(51) 



Using the relations (IA.5I) . the coefficient c can be ex- 
pressed in terms of the Debye mass as 



24tt 2 



(N 2 c ~\)Lm 2 D . 



(52) 



Furthermore, in equilibrium, a = cT. 

We note that in spite of our neglect of transverse chro- 
modynamic fields, the collision terms for the isotropic 
plasma derived here are very similar to those derived in 

til mil Hi 11^3. 



C. Equilibration of an isotropic plasma 

As an application of the Fokker-Planck equations (|38|) 
we discuss the problem of plasma equilibration. In this 
section we limit our considerations to quarks, as the anal- 
ysis for antiquarks and gluons is very similar. We con- 
sider the system which is homogenous and mostly equi- 
librated but a small fraction (A <C 1) of the particles, 
denoted by Sn(t,p), is out of equilibrium. One asks on 
what time scale the system reaches the equilibrium. The 
distribution function is assumed to be of the form 



n(t, p) = (1 - A) n cq (p) + A 5n(t, p) 



(53) 



In the course of equilibration nit, p) tends to n cq (p). 
Since the particle number is conserved within the trans- 
port theory approach developed here, 8n(t, p) is not re- 
duced to zero in the equilibration process but it tends to 
n eq (p). 

We define the rate of equilibration T through the rela- 
tion 



I- 1 *- 



(54) 



We note that T is either positive, when Sn grows going 
to n cq , and it is negative, when Sn decreases going to 
n cq . Using the Fokker-Planck equation (j38]), the defini- 
tion (|54"|) gives 



r = 



1 



y i p X ij V 3 p + %Y l )Sn 



(55) 



Since the fraction of particles with non-equilibrium dis- 
tribution is assumed to be small, the coefficients a, 6, c 



from the formulas (|44l 1451) are given by the equilibrium 
function n cq ~ e~ Ep / T . Using the approximate expres- 
sion of a (|50| with 6 = and c = a/T, Eq. (|55|) is rewrit- 
ten as 



TE r 



5n . 



(56) 



The equilibration rate obviously depends of the form 
of Sn. Here we consider the case where the small fraction 
of partons has an equilibrium distribution of temperature 
To which differs from the temperature T of the bulk of the 
partons. Thus, Sn ~ e _Bp ' T °. Then, the equilibration 
rate (|56|) equals 



T = a 



T-T 
T 2 TE P 



(E P - 2T ) 



(57) 



For T — To, the whole system is in equilibrium and, as 
expected, V = 0. When T > T , the distribution e~ Ep/T ° 
is steeper than e~ Ep ^ T . Equation (|57|) tells us that Sn de- 
creases for E p < 2Tq and grows for E p > 2To during the 
equilibration process. When T < T), we have the oppo- 
site situation. In both cases, the slope of the distribution 
function Sn tends to the slope of n cq . With the coefficient 
a given by (fST))) . the formula ([57]) quantitatively predicts 
how fast the equilibrium is approached. 



VI. TWO-STREAM SYSTEM 

The two-stream configuration provides an interesting 
case of an unstable plasma. The correlation function of 
longitudinal chromoelectric fields, which is needed to de- 
rive the transport equations, was computed in (43[. Un- 
fortunately the correlation function for transverse fields 
is not known. This limits our considerations to longitu- 
dinal fields. 

The distribution function of the two-stream system is 
chosen as 



f(p) = (27rfn \S^ (p - q) + <5 (3) (p + q) 



(58) 



where n is the effective parton density in a single stream. 
The distribution function (|58p should be treated as an 
idealization of the two-peak distribution where the par- 
ticles have momenta close to q or — q. 

To compute £r,(w, k) we first perform an integration 
by parts in (|A.1|) and then substitute the distribution 
function (|5"5|) into the resulting formula. We obtain 



9 



£i(w,k) = 1 - 



Jk-u) 

k 2 



2 r 



(w - k • u) 2 (uj + k • u) 2 



(59) 



_ (w - w+ (k)) (u> + lo + (k)) (u - u>- (k)) (w + w_ (k)) 
~ ( W 2 -(k-u) 2 ) 2 ' 

where u = q/-E q is the stream velocity, /i 2 = g 2 n/2Eq and ±w±(k) are the four roots of the dispersion equation 
Ei(ct),k) = which are explicitly given by 



c4(k) 



k 2 



k 2 (k • u) 2 + /i 2 (k 2 - (k • u) 2 ) ± ii J(k 2 - (k • u) 2 ) (4k 2 (k • u) 2 + ^ 2 (k 2 - (k • u) 2 ) 



(60) 



One can show that < w+(k) <E R for any k, while w_(k) is imaginary for k-u 7^ and k 2 (k-u) 2 < 2/i 2 (k 2 — (k-u) 2 ). 
uj- represents the well-known two-stream electrostatic instability generated by a mechanism analogous to the Landau 
damping. For k 2 (k • u) 2 > 2^i 2 (k 2 — (k • u) 2 ), the cj_ mode is stable: < w_(k) G E. 

The terms like (E(t, r) • V p SQ(t, r, p)), which enter the transport equations (|T2"j) are given by Eqs. (TT^l) . As for the 
isotropic plasma one needs to specify the the correlation functions (w, k)i?£(a/, k')), (E(t, r)^Qo( r ', p'))> etc - The 
correlation function of the longitudinal fields k)^(w',k')) was found in 



<£>,k)i^(a/,k')> = -. 9 2 <5 



2XQ , (27r) 3 5( 3 )(k + k')^F 



k 2 k 2 

• 2 - (k • u) 2 



(k-u)(k'.u) 



(61) 



to 



(uj - w_(k)) (w + w_(k)) (w - w+(k)) (u> + w+(k)) 



(k' ■ u) s 



(w' - w_(k')) (w' + w-(k')) (w' - w+(k')) (w' + w+(k')) 



We are particularly interested in the contributions of the unstable modes to the correlation function. For this reason 
we consider the domain of wave vectors obeying k • u ^ and k 2 (k • u) 2 < 2/s 2 (k 2 — (k • u) 2 ) when uj- (k) is imaginary 
and the mode is unstable. We write w_(k) = 17k with < 7k G E. The contribution coming from the modes ±w_(k) 
then equals [43| 



(K(t,r)^(t',r')) U nst a ble = ^S ab n 



d 3 k e lk ( r 



(7^ + (k-u) 2 ) : 



(2tt) 3 k 4 (uj 2 



(62) 



( 7 £ + (k • u) 2 ) cosh ( 7k (< + t')) + ( 7 £ - (k • u) 2 ) cosh ( 7k (< - t')) 



As Eq. ((HI]) shows, the contribution of the unstable 
modes to the field-field correlation function is space 
translation invariant - it depends only on the difference 
(r — r'). If the initial plasma is on average homogeneous, 
it remains so over the course of its evolution. The time 
dependence of the correlation function (|62[) . however, is 
very different from the spatial dependence. The elec- 
tric field grows exponentially and so does the correlation 
function, both in (t + t') and (t — t'). The fluctuation 



spectrum also evolves in time as the growth rate of the 
unstable modes is wave- vector dependent. After a suffi- 
ciently long time the fluctuation spectrum will be domi- 
nated by the fastest growing modes. 

The correlation function (E(t, r)SQo(r — vt, p)) is, as 
previously, given by Eqs. (|29[) . Since the dielectric func- 
tion (|59"|) is real, the correlation functions (E(t, r)5Qo(r — 
vt,p)>, <E(t,r)*Q (r-vt,p)> and (E(t,r)SG (r-vt,p)> 
all vanish. Therefore, 



Tr(E(t,r).V p <5Q(t,r,p)) = gf dt' V* r)i^(i',r - v(t -i'))>V>(p) 

Jo 
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9 {Nt-l)n I dt'% 



d 3 k e ikv (*- t ') WW ( 7 £ + (k-u) : 



i y " c ~ J ~J a ~~ ' P J (2tt) 3 k 4 (w 2 -w 2 ) 2 
x + (k • u) 2 ) cosh ( 7k (t + 0) + (7k - (k • u) 2 ) cosh ( 7k (t - t'))\ Vjn(p) . 
Performing the integration over t' and keeping only the real part, one finds 



(63) 



Tr(E(i,r)-V p 5Q(i,r,p)> = ^ (iV c 2 - 1) n V* 



d 3 fc 



(7^ + (k-u) 



(27T)3 k 4 ( W 2 -C 2 )2 7 2( 7 2 + (kv) 

7 £ sinh(2 7k t) + (k • u) 2 ( 7k sinh(2 7k i) 



(64) 



+ (k • v) sin ((k • v)t) cosh( 7k i) — 7k cos ((k • v)t) sinh( 7k i) 
Neglecting the oscillating terms, we finally get 

d 3 k WW 



V>( P ) 



Tr(E(t,r)-V p <5Q(i,r,p)) 



9 -{N 2 c -l)nV\ 



( 7 2 + (k-u) 2 ) 3 



(27T)3 k 4 (c 2 - W 2 ) 2 7k ( 7 2 + (kv)2) 



sinh(2 7k t) V>(p) . (65) 



In an analogous way one can obtain explicit expressions for (E(t, r) • V p SQ(t, r, p)) and (E(t, r) • V„6G(t, r, p)). We 
do not present these here, because they do not provide any new insight. 

Since we explicitly integrated over the distribution function (|58l) in deriving these results, we only give the transport 
equations (fT2"j) for the two-stream system in the Fokker-Planck form: 



id- v;x''./.v)v;,)^^;P^=«, 



(©-VjX«(t,v)V£K(t,r,p) = 



where 



WW ( 7k + (k • u) : 



g 4 N 2 — 1 f d 3 k 
XV{t ' V) = 1^N~ n J (2^kV 2 + 7 2 ) 2 7k ( 7 2 + (kv)2 



sinh(2 7k i) , 



(66) 



(67) 



and X*j(t,v) = 2N 2 X^(t,p)/(N 2 - 1). 



To get an idea how the two-stream system evolves ac- 
cording to the Fokker-Planck equations (|6"6"|) . we take into 
account in the integral (|67[) only those wave vectors which 
are parallel to the stream velocity u, with the latter being 
chosen along the axis x. The only non- vanishing compo- 
nent of X 13 (t, v) is then X xx (t, v). Neglecting the depen- 
dence of X xx (t, v) on p and assuming that the system is 
homogenous, the Fokker-Planck equation (|66|) for quarks 
becomes a one-dimensional diffusion equation 



dn(t,p) D ^d 2 n(t lP ) 



dt 



dp 2 x 



(68) 



with the diffusion coefficient D(t) = X xx (t) depending 
on time approximately as 



(69) 



D(t) = de 

where d and 7 are constants. 

If the distribution function is initially of the form 

n(* = 0,p)=27rn*(p 1B -g) J (70) 



where n is independent of p x , the solution of the diffusion 
equation (|6"8"]l is found as 



n(t,p) 



27T 7 



d(e 2 ^ - 1) 



exp 



l{Px - qf 



2d{e 2 ^ - 1) 



(71) 



The distribution function (|7ip is normalized in such a 
way that 



dp x 
2vr 



n(t, p) = ft . 



According to the solution (I7T1) . the electric field growing 
due to the electrostatic instability rapidly washes out the 
peak-like structures of the two-stream distribution func- 
tion (|S5J). It should be understood, however, that the 
solution (fTTj) is valid only for time intervals which are 
sufficiently short that the distribution function used to 
compute the coefficient X t3 (t,v) is not much different 
from the function (|5"5|) . Nevertheless, the solution (fTTj) 
shows how the equilibration process commences. 
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VII. SUMMARY AND OUTLOOK 

Wc have developed here the quasi-linear transport the- 
ory of a weakly coupled quark-gluon plasma. Our main 
motivation was to study the equilibration of plasmas that 
are initially unstable. The field fluctuation spectrum, 
which is found within the linear response approach, de- 
termines the evolution of the regular distribution func- 
tions. More specifically, the fluctuations of chromody- 
namic fields provide collision terms to the transport equa- 
tions of the regular distribution functions. We have 
limited our considerations to longitudinal chromoelectric 
fields, as then the field correlation functions are known 
for both the isotropic and two-stream systems. The col- 
lision terms were found in either the Balescu-Lenard or 
Fokker-Planck form. In the case of an isotropic plasma we 
showed how the system equilibrates when a small fraction 
of particles has a different temperature than the bulk. 

The case of the two-stream system is more interest- 
ing. The Fokker-Planck equation could be approximately 
written as an equation of diffusion in momentum space. 
The diffusion coefficient, which is given by the chromo- 
electric fields for the two-stream instability, exponentially 
grows in time. We found the exact solution of the diffu- 
sion equation, which showed that the peak-like structures 
in the parton momentum distribution dissolve rapidly. 

In nonrelativistic plasmas it is often a well justified ap- 
proximation to keep only longitudinal electric fields and 
to neglect magnetic and transverse electric fields |4l[ . 
In the case of ultrarelativistic plasmas, this is no longer 
true. If initially the fields are purely longitudinal, the 
transverse fields are automatically generated, and they 
are dynamically important. Therefore, the ultrarelativis- 
tic plasma considered here, where the transverse fields 
are neglected, should be rather treated as a toy model 
which we have studied mostly for the sake of analyti- 
cal tractability. With this simplified example we have 
been able to elucidate some general features of the prob- 
lem. Physically better motivated situations will require 
substantial numerical work, which is less conducive to 
general insights. 

The considerations presented here clearly demonstrate 
the usefulness of the quasi-linear transport theory for 
the study of equilibration processes of quark-gluon plas- 
mas. As mentioned in the Introduction, numerical stud- 
ies indicate that the unstable chromomagnetic plasma 
modes play an important role at the early stage of the 
quark-gluon plasma produced in relativistic heavy-ion 
collisions. Therefore, it would be of considerable inter- 
est to compute the correlation functions of transverse 
fields in arbitrary anisotropic plasmas in order to de- 
rive the relevant transport equations. As explained in 
[43| . there is no conceptual difficulty in such a compu- 
tation, but one has to invert the matrix E IJ (u;,k) = 
— k 2 ^ + k l ti + uj 2 e r} (w,k). This is easily done for 
isotropic plasmas but for anisotropic plasmas one obtains 



a rather complex expression which is very cumbersome 
for further analytic calculations [46]. Except for some 
special cases, numerical methods seem to be unavoid- 
able. Such computational studies are beyond the scope 
of the present work but progress in this direction will be 
hopefully reported soon. 
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Appendix 

We discuss here the longitudinal chromodielectric per- 
meability el(w, k) which is known to be 



<r / d 3 P k-v p /( P ) 



2k 2 J (2tt) 3 w-k-v + ifU 
Applying the identity 

„ 1 1 - n+ = V-^m5{x) 



(A.l) 



to Eq. (IA.1I) . one immediately finds Syex,(o;,k) 

3e L (c, k) = - J ^ 27n5( W - k • v) k • V p /(p) . 

(A.2) 

If the plasma is isotropic V p /(p) can be expressed as 

(A.3) 



Y7 ft ^ d f(^ 
V p /(p) = ^T^V 



dE„ 



And if the partons are additionally masslees, the integral 
in (|A.ip factorizes into the angular integral and the inte- 
gral over p = |p|. Then, one finds the real and imaginary 
parts of the longitudinal chromodielectric permeability 
£l(w, k) as 



1 — n 

2 k 



w-|k| 



sMo,,k) = |e(k 2 - w 2 )^ 



where the Debye mass tud is 
a 2 

,2 _ y 



m D = -- 



An 



dpp 



2 d/(p) 







dp 



(A.4) 



(A.5) 
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